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Abstract 

We make a thermal analysis of the light scalar susceptibilities using SU(3)-Chiral Perturbation 
{SJ Theory to one loop taking into account the QCD source of isospin breaking (IB), i.e corrections 

j> coming from m u ^ m&. The value of the connected scalar susceptibility in the infrared regime, the 

one relevant when approaching chiral symmetry restoration, and below the critical temperature 
is found to be entirely dominated by the 7r° — rj mixing, which leads to model-independent O(e ) 
£T) corrections, where e ~ ma — m u , in the combination \ uu — Xud of flavour breaking susceptibilities. 

We also present preliminary results for the corrections to the real part of the pion self-energy at 

, , nexst-to-leading order in SU(2)-Chiral Perturbation Theory taking into account electromagnetic 

i-^ interaction. The zero and finite temperature results for the charged and neutral pions are given in 

^ H terms of the three-momentum of the external pion; and their difference is calculated to this order 

^-i stressing the fact that, at low and moderate temperature, the mass splitting M n ± — M n o grows 

with temperature for, at least, non-zero charged pion mass running inside the loops. 



> 1 Introduction 

The low-energy sector of QCD has been successfully described over recent years within the chiral 
*0 lagrangian framework. Chiral Perturbation Theory (ChPT) is based on the spontaneously breaking 
of chiral symmetry SUv(Nf) x SUA(Nf) — > SUv(Nf) with Nf = 2,3 light flavors and provides a 
^~~] consistent, systematic and model-independent scheme to calculate low-energy observables [HEIE]. The 
t-h effective ChPT lagrangian is constructed as an expansion of the form C = C p 2 + £ p 4 + . . . where 
>• p denotes a meson energy scale compared to the chiral scale A x ~ 1 GeV. The formalism can also 
^ e ex t en ded to finite temperature T, in order to describe meson gases and their evolution towards 
5-h chiral symmetry restoration for T below the critical temperature T c [U [5] , where T c ~ 180-200 MeV 
from lattice simulations P, E]- The use of ChPT in this context is important in order to provide 
model-independent results for the evolution of the different observables with T, supporting the original 
predictions for chiral restoration [9], also confirmed by lattice simulations, which are consistent with a 
crossover-like transition for Nf = 3 (2+1 flavours in the physical case), a second-order one for Nf = 2 
in the 0(4) universality class and first order in the degenerate case of three equal flavours. 

The invariance under SUy(2) vector group is the isospin symmetry, which is a very good approxima- 
tion to Nature. However, there are several proccesses where isospin breaking corrections are phenomeno- 
logically relevant, for example sum rules for quark condensates [3j, meson masses [TU] or pion scattering 
[Tl~1, [T2"]. There are two possible sources of isospin breaking: the QCD — m u light quark mass differ- 
ence and electromagnetic interactions. Both can be accommodated within the ChPT framework. From 
the first source we expect corrections of order (rrid — m u )/m s , encoded in the quark mass matrix, which 
generates also a 7r°r] mixing term in the SU(3) lagrangian [3]. On the other hand, the electromagnetic 
interactions are included in the ChPT effective lagrangian via the external source method and give rise 
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to new terms [TUl [HI [121 EH] of order C e 2, C e 2 p 2 and so on, being e the electric charge. It is possible to 
accomodate these terms into the ChPT power counting scheme by considering formally e 2 = 0(p 2 /F 2 ), 
with F the pion decay constant in the chiral limit. 

The purpose of this paper is to calculate the leading thermal contributions to the connected and 
disconnected scalar susceptibilities taking into account isospin breaking, and to show our preliminary 
results about the thermal evolution of the masses of the charged and neutral pions. 



2 Light scalar susceptibilities and the role of the tt° — 77 mixing 

Light quark condensates, (uu) and (dd) , at zero and finite temperature have been calculated to one loop 
in [TJ] and [15] . respectively, within the framework of SU(3)-ChPT taking into account both sources of 
IB. The main feature we want to stress is that there is a it — rj mixing term appearing through the 
tree- level mixing angle e defined by tan2e = ^ m ' i ~ m " . 

Different light quark masses allow to consider three independent light scalar susceptibilities defined 

as 

X* = = 9^. ^Z(m u * m d ). (1) 

From now on in this section, we will neglect the electromagnetic corrections because they are small 
and not relevant for our present discussion, so we will put e = 0. Then, to leading order in the mixing 
angle, the contribution of the 7T°—r] mixing in the quark condensate sum is of order e 2 whereas for the dif- 
ference it goes like e. The thermal functions ^(T, Mj), i = 7r°, 77; defined as^(T) = 4n l F2 / °° dp^- $e p _ v 
with E 2 = p 2 + M 2 and (3 = T _1 ; are suppressed by those coefficients and the quark condensates do 
not receive important corrections. The important point is that differentiating with respect to a light 
quark mass is essentially the same as differentiating with respect to e ~ nid ^ l rnu , so the suppression of 
the thermal functions is smaller in the case of the susceptibilities than in the quark condensate. 

Because of the linearity in e in (uu — dd) for a small mixing angle, the combinations Xuu ~ Xud and 
Xdd — Xdu receive an 0(1) IB correction due to 7r° — 7] mixing, which would not be found if m u = rrid 
is taken from the beginning. The analysis of the e-dependence of (uu — dd) shows that, up to 0(e), 
Xuu — Xdd, so combinations like Xuu — Xdd, which also vanish with m u = md, are less sensitive to IB. 

One can also relate these flavour breaking susceptibilities with the connected and disconnected ones 
US], often used in lattice analysis PU US]: Xdis = Xud, and Xcon = \(Xuu + Xdd - ^Xud)- From the 
previous analysis, we get Xcon - Xuu - Xud- 

Therefore, our model-independent analysis including IB effects provides the leading nonzero contri- 
bution for the connected susceptibility which arises partially from 7r° — 77 mixing. This is particularly 
interesting for the lattice, where artifacts such as taste-breaking mask the behaviour of Xcon with the 
quark mass and T when approaching the continuum limit [H]. In fact, our ChPT approach is useful 
to explore the chiral limit (m u ^ — > 0) or infrared (IR) regime, which gives a qualitative picture of the 
behaviour near chiral symmetry restoration. In this regime M n < T <C Mk, and therefore we can 
neglect thermal heavy particles, which are exponentially suppresed. 

The leading order results for the connected and disconnected susceptibilities at zero temperature 
are the following 

xit(T = 0) = 8B 2 [2L r s (f,) + H*M]) - -±- 1 + log -f - -JL log -± + 0(e 2 ) (2) 
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X l«(T = 0) = 325 2 L^) - ^1 + log^J + 2^ (5\og^ - 1 ) +0(c*). (3) 



2 




200 



Figure 1: Connected IR susceptibility normalized to B%, for fixed tree level eta mass and m s 



80 MeV. 



where B is the parameter which relates masses and quark condensates at tree level via the Gell 
Mann-Oakes-Renner formula, and L 6 , L 8 , H 2 are low energy constants. 

The log term in equation ^ is the dominant at T = and can be found in |l6j . but the connected 
IR susceptibility ^ is not zero at T = 0, because it receives contributions of order 0(1) in the mixing 
angle. 

If we consider the pion gas in a thermal bath, then expressions — (J3j) are modified according to 
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Note that, as we have already mentioned, the eta mass term in equation (|4]) and in the subleading 
corrections in the mixing angle comes from the e- analysis and the IR expansion of the g\ (M n ), and does 
not have anything to do with thermal etas. 

Figure [I] and [2] show, respectively, the connected susceptibility Q for fixed tree level eta mass 
(proportional to y/Bo m s in the IR regime), and the disconnected one (pi) for several values of the light 
quark mass ratio m/m s , and also with fixed tree level eta mass. The leading scaling with T and the 
light quark mass in this regime for the disconnected piece goes like i.e the same scaling calculated 
in (THl [T7j; whereas the connected susceptibility grows quadratic in T over a mass scale much greater 
than the SU(2) Goldstone boson's one. Therefore, in the continuum limit, we only expect x<n.s to peak 
near the transition, as the m — > + limit in figure clearly shows. 



3 Thermal charged and neutral pion masses in SU(2)-ChPT 

at 0{p A ) 

If we consider virtual photons in the calculation of the real part of the self-energy in the mass shell, there 
are four relevant diagrams that correct the masses at order 0(p A ): there are pion tadpoles, diagram(a) 
in Figure [3j and the tree level NLO diagram needed for renormalization, (b), where both charged 
and neutral pions participate; and diagrams with virtual photons, (c) and (d), which only modify the 
charged pion mass. 

The photon-tadpole diagram (c) is proportional to the photon mass and therefore vanishes at zero 
temperature, while pion tadpoles (a) and the photon-exchange diagram (d) are finite and chiral scale- 
independent once regularized and combined with diagram (b). The LO corrections to the masses of the 



3 




50 100 150 200 

T(MeV) 



Figure 2: Disconnected IR susceptibility normalized to Bq, for several light quark mass ratios and fixed 
tree level eta mass (m/m s = 0, 05 is the physical case). 




Figure 3: Diagrams for the pion self-energy at order 0(p 4 ). (a), (b) represent pion tadpole contributions 
and the tree level NLO diagram wich renormalizes the loops, respectively (c) (photon tadpole) and (d) 
(one-photon exchange) only correct the charged pion self-energy. 



SU(2) NGB are calculated as M 2 = M 2 + E(M 2 ), where M is the respective tree level mass. At zero 
temperature the neutral and charged pion masses are [11] 



M 2 = M 2 1 + 2fi ± -n + eXo + 21 



M 2 \ B 2 4 

) - 2—l 7 (m d - m u ) 2 - -Be 2 k 7 (m d - m u ) (6) 



where /i ±i0 = 3 ^ 2 f 2 log , and /C^o = -f [k[ + k r 2 - ^(2^ - k\) - k\ - k r 6 - \k r 7 ] , being h elec- 

tromagnetic low energy constants, Z the parameter that corrects the leading order charged pion mass; 
and 
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1 + — + hq + e 2 /C^± + 21 
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(7) 



with the definitions 
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and 
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~9~ 



2Z (k[ + ¥ 2 ) - -k 13 - k u 



There is a factor 1/2 in the coefficient of the kj for both masses wich does not appear in [TT], and 
that was also noted by [i~9] . 
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To this point the total mass difference between charged and neutral pions becomes 



B 2 

M 2 ± - M 2 = 2M 2 (/i - M±) + ^l r r (m d - m u f + 2e 2 F 2 



+e 2 )Ci} + M> 2 



f 

47T 



-2e 2 F 2 (3 + 4Z)fi±, (8) 



which has pure strong, pure EM and mixed EM-strong contributions. 

If the pions are immersed in a thermal bath, there appear new contributions with no UV divergences, 
since these only appear in the zero temperature part and they have been already renormalized. The 
pion tadpoles, charged or not, give rise to g%(M,T) thermal functions. For the neutral pion mass we 
get 



M 2 = M 2 o {T = 0) 



l + iL(M^,T)-i ?1 (^„,r) 



(9) 



where is worth noting that the neutral pion mass decreases with T, contrary to what happens if we do 
not consider electromagnetic effects, as can be seen if we put M^ ± = M 2 in the last expression. 

As for the charged pion we can separate the contributions coming from pion tadpoles, M 2 tadpoles ; and 
the two different contributions from the virtual photon diagrams: one coming from the photon tadpole 
which is not zero at finite temperature and gives a typical thermal screening contribution, Mp h tadpole ; 
and the other due to the one-photon exchange, M| h Exchange . 

M 2 

For the first to ones we get M 2 Tadpoles = ^ -4Ze 2 9l (M 2 ± , T) and M 2 h Tadpole = \e 2 T 2 . The latter 
has the typical form of a Debye or screening mass of the electric field in a thermal bath |20] , which 
always grows with T. 

The one-photon exchange diagram is more complex and its contributions to the real part depend, 
in general, on the three-momentum of the external pion, which is a direct consequence of the Lorentz 
symmetry breaking in the thermal bath. The Matsubara sums can be performed in the standard way, 
before performing the analytic continuation to external continuous frequencies 



., . . d 3 k n(u) (2q- k) 2 \ k()=UJ f d 3 k n(u) (2q - k) 2 

tie (^p n . Exchange J ~ e \ / \ 3 ~7, 7 ^ ^ + 



fco= 



+ 



(2vr) 3 2tu (qo-u) 2 -u' 2 J (2tt) 3 2u (q + u) 2 - uo' 2 
d 3 k n(w') (2g - /c) 2 | fco =g _a;' , f d 3 k n(u') (2q - k) 2 \ ko=qo+u) < 



(2tt) 3 2co' (g - oj') 2 - cu 2 J (2vr) 3 2cu' (g + to') 2 - u 



with bj 2 = k 2 the photon energy squared inside the loop and u' 2 = (q — k) 2 + M 2 ± the pion energy 
squared, also inside the loop; being k the virtual photon four-momentum. The above three-momentum 
integrals can be written as 

e 2 f°° r „ s M 2 ± JE 2 - M 2 ± cos - k(E 2 sin 2 + M 2 ± cos 2 0) , s 

Re(E Ph . Exchange = — / dk tyn (k sin0^^ g . 2 \ 1 f ^ ^ 10 

^ Jo Jo E sm r + cos 4> 

e 2 f°° r n(cu') k 2 (M 2 ± cos 2 + E 2 sin 2 0) + 2M 2 ± (E 2 - kJE 2 - M 2 ± cos 0) 
' dk I do sin ^ 



4vr 2 Jo Jo u' T E 2 sin + M^ ± cos 2 

It is clear now that the charged pion real part of the self-energy depends on the energy of the external 
pion E 2 = |g| 2 + M 2 ± , and, therefore, on the external three- momentum. 

The final result will take the form M 2 ± = M 2 ±(T = 0) + M 2 Tadpolcs + M 2 h Tadpole +Re (S 

Ph. Exchange ) ■ 

In Figure [4] we have plotted our preliminary results for the neutral and charged pion real masses as 
a function of the temperature of the thermal bath, taking physical values for the pion masses and 
calculating Mp hExdl e in the static limit, i.e, for E = M n ±. Once we have the thermal and isospin 
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Figure 4: Preliminary results for the charged (dashed line), neutral (solid line) masses and its difference 
(dot-dashed line) at LO in the static limit. 

breaking corrections to the masses separately for any value of the external momentum, we can calculate 
it in the limit where temperatures are (i) much greater than the masses and the external momenta 
(which means that we have to set the masses inside the loops to zero), and (ii) sizeable to the momenta 
running inside the loops, T ~ k 3> m, q . With this assumptions, we are lead to the HTL result given 

in [21], M%± — M^ = M%± ^1 — =^rj + |e 2 T 2 , which serves us as a consistency check. 

Moreover, our low temperature analysis allows to assume a slightly different chiral limit, in the sense 
that we can still assume m u = = but considering e ^ even inside the loops. In Figure [5] we 
have plotted our preliminary calculation both in this latter limit, and also considering m u = rrid ^ 0, 
e 7^ 0; to be compared with those appearing in the HTL result [21] where the contributions from the 
screening-like term, always increasing with T and inherent to the thermal bath, are responsible for 
the final growth of the mass difference, contrary to what one would expect naively from the sum rule 
relating the axial and vector spectral functions [20] applying chiral symmetry restoration arguments. 
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Figure 5: Different results for the charged-neutral pion mass difference, (a, solid line) corresponds to 
our preliminary results in the chiral limit keeping corrections e ^ for the tree level charged pion mass 
inside the loops, (b, dot-dashed line) corresponds to the same preliminary calculation with m u = ^ 
and e^0 also inside the loops; and the full dashed line is the result given in [2"T| . 
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